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Periodicity and relaxation are investigated for the trajectories of the states in 
cylindrical linear cellular automata. The time evolutions are described with 
matrices. The eigenvalue analysis is applied to obtain the maximum values of 
period and relaxation. The translational invariance suppresses the maximum 
period and relaxation. 
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I. INTRODUCTION 



Cellular automata (CA) are one of the simplest mathematical models for nonlinear dy- 
namics to produce complex patterns of behavior. They had been originally introduced by 
von Neumann [[]J to investigate some artificial life. Wolfram had reintroduced CA as a model 
to investigate complexity and randomness || . He investigated many fundamental features of 
them Since then many authors have made efforts to clarify the properties of CA and 

applied to natural systems including biological systems, complex fluids, chemical reactions, 
astronomical systems and so on ||. 

Computational simplicity is the advantages of using CA to simulate complex behaviors. 
Discrete states are defined on each lattice point and evolve with discrete time steps in CA. 
Moreover nonlinear properties are expected to be emphasized by discretization. A simple 
modeling, however, does not mean simple properties. Even though the great efforts, the 
fundamental properties of CA are still unclear. Studies on these simple systems seems to be 
a key for understanding discretized simulated systems. 

One-dimensional CA are described by the discrete time evolution of site af 

Oi(t + 1) = F[a,i- r (t), Oi_ r+1 (t), . . . , Oi(t), . . . , a i+r (t)), (1.1) 

where a, takes k discrete values over Z^. The simplest models, elementary cellular automata, 
consist of sites with two internal states over Z2 (k = 2) interacting with the nearest neighbor 
sites (r = 1). Wolfram introduced a naming scheme for these models and classified the 
behavior of CA into four classes @|§. 

Most authors have worked on CA within the scope of the infinite number of sites. A few 
works have concerned the effects of finiteness. The orbits of cylindrical linear CA, which are 
CA with periodic boundaries, are analyzed on the basis of characteristic polynomials, which 
directly describe the states of linear CA [pj,fT|— [T0[1 - Stevens, Rosensweig and Cerkanowicz had 
investigated rule 90 CA, which will be described in ||, with Dirichlet boundary conditions. 
They analyzed the eigenvalue polynomials of the matrices which give the time evolution of 
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the system [pjj . In our previous papers fll"2"l , |13"|l (referred as papers I and II), we had also 



investigated the periodic orbits of finite linear CA (rule 90 and 150) with Dirichlet boundary 
conditions by analyzing the eigenvalue equations. In the present paper the method is applied 
to cylindrical linear CA. The proof of the classification of the orbits for Dirichlet boundary 
cases is given in Appendix [A]. 

II. THE MODELS AND ORBITS OF STATES 

There are two examples of linear cellular automata in the elementary ones (k = 2, r = 1). 
They are called as rule 90 and rule 150 following the Wolfram's naming scheme. For rule 
90 CA, the time evolution of the i-th site a« £ {0, l}(i = 1, . . . , N) is described as a sum 
modulo 2 of the nearest-neighbor sites: 

a,i(t + 1) = ai_i(t) + a i+ i(t) mod 2. (2.1) 

The time evolution of rule 150 CA is given as a sum modulo 2 of the nearest-neighbor sites 
and itself: 

di(t + 1) = Ctj_i(t) + CLi(t) + Oj+l (t) mod 2. (2.2) 

We use the periodic boundary conditions: a = a^, a/v+i = a%. These models belong to the 
third class which show the chaotic behavior in the Wolfram's classification 

These models are linear because the rules are additive and the time evolutions are also 
expressed by the matrices 

A(t + 1) = UA(t), (2.3) 

where A(t) = t (ai(t), ci2(t), . . . ,aisf{t)) describes the state, N bits binary number, at t. The 

components of the transfer matrices U are given by 

f 1 j = i±l 

1 i = 1, j = N 

(rule-90), (2.4) 

1 i = N, j = 1 
* otherwise 
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Uij = < 



' 1 


j = i ± 1 




1 


j = i 




1 


i = 1, j = 


iV 


1 


i = N, j = 


-- 1 


- 


otherwise 





(2.5) 



The trajectories of the states can be classified into three cases with the properties of the 
transfer matrices as mentioned in papers I and II (See Fig. |]). 

U Un = I, (2.6) 

JjU N +TT N _ JJ^N (2 7) 

U fN = 0, (2.8) 

where I denotes a unit matrix. The first case, Eq. (|2.6|) , corresponds to perfectly periodic 
motions with the maximum period Un- Namely all states belong to perfectly periodic orbits 
whose period does not exceed II at. The second case, Eq. ( |2.7| ), shows periodic motions with 
relaxation. A relaxation path whose time step is not greater than ir^ leads to a periodic 
orbit with the maximum period Ily- The third, Eq. ( |2.8j ), is a special case of the second 
one. All states are drawn to a null state on the contrary to that the current CA belong to 
the third class of Wolfram's classification of motions (chaotic motions). 



III. EIGENVALUE ANALYSIS 

The eigenvalue polynomials for A-site linear cellular automata are defined by 

D$\\) = \U + XI\, (3.1) 

where the index R specifies the rule number and the boundary condition, for example, 90D 
for rule 90 CA with Dirichlet boundaries and 150P for rule 150 with periodic boundaries. 
Since each site a« takes binary values, the eigenvalue polynomials are over Z2, namely each 
coefficient is or 1. The eigenvalues A are not usual numbers but over the Galois field 
GF(2 Af ), a finite field with 2 N elements. 
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The eigenvalue polynomials enable us to find the maximum period Iljv and the maximum 
relaxation as discussed in Ref. [IT], papers I and II. For nilpotent cases, Z)jy (A) = X N , 



all states are drawn into a null state within iV steps or less. If the eigenvalue polynomial 
has a constant term, D N (0) = 1, the eigenvalue equation D N (A) = can be reduced to 
a simple form as A Pjv + 1 = by multiplying some powers of A and repeatedly substituting 



the eigenvalue equation ||14|| . Then the minimum value of Pjy corresponds to the maximum 
period Hat. In other words the maximum period Un is the minimum integer m satisfying 
D N (X)\(X m — 1), where the notation f\g denotes / divides g. The remaining cases are 
mixture of above two cases as D N (A) = X PN x d(X) , where the polynomial d( A) has a constant 
term (D$\o) = and d(0) = 1). One can evaluate the maximum period with the above 
procedure from d(X). The maximum period is the minimum integer m satisfying <i(A)|(A m — 
1). The relaxation length to periodic orbits is p^. Actual periods and relaxation paths are 
given by divisors of the maximum values respectively. They depend on the symmetry of the 
initial states as shown in papers I and II. 

The eigenvalue polynomials for rule 90 cylindrical CA obey the relation 

< P) (A) = A<^(A). (3.2) 

Orbits of rule 90 cylindrical CA can be classified as 

u ii N +i = u (AT i s dd), (3.3) 
ljn N +ir N = jj7r N ( N ig even except n = 2 n ), (3.4) 
U N = (N = 2 n ), (3.5) 

by Eq. (|3.2|) and the classification of orbits for rule 90 finite CA with Dirichlet boundaries 
(see Appendix 0). There are no perfect periodic orbits by comparing with the Dirichlet 
boundary cases. 

The maximum values of the period and relaxation are derived from the polynomials 

L(JV-l)/2j 

Df P) (A)= J2 (C?- l mod2)X N - 2 i, (3.6) 

3=0 



which are the solutions of Eq. ( |3.2|) , where [^J denotes the largest integer not exceeding x 
and 

C?={-iy{^~ j y (3.7) 

Some examples of the eigenvalue polynomials are shown in Table |I[ Equation ( |3.2|) leads that 
the maximum period for rule 90 cylindrical CA with N sites are equal to that for A" — 1-site 
case with Dirichlet boundaries, ITv(rule 90P) = Il A r_ 1 (rule 90D). 

The eigenvalue polynomials for rule 150 cylindrical CA can be also written with those 
with Dirichlet boundaries 

^ 50P) (A) = (l + A)DrT ) (A). (3.8) 

Orbits can be classified as 



U n " = I 


(N^3n), 


(3.9) 


U N = I 


(N = 2 n ), 


(3.10) 


JjU N +TT N _ JJKN 


(N = 3n). 


(3.11) 



by virtue of the classification of orbits for rule 150 CA with Dirichlet boundaries. There are 
no characteristic difference from the Dirichlet boundary cases. No simple relations on the 
maximum periods between the periodic and Dirichlet boundary cases are found. 
The maximum periods and relaxations are derived from the polynomials 



'L(JV-1)/2J \ N (\_{N-k)/2\ /M 0n\ 

^ &UFj (A)=l E ^- 1 mod2 +E E Cf~ x ( J mod 2 I A fe . (3.12) 



h 

3=0 J k=l \ j=0 \ ft 

Some examples of the eigenvalue polynomials are shown in Table [TT]. 

The actual period and relaxation of the system may be suppressed by the translational 
symmetry which the periodic boundary conditions ensure. Direct matrix multiplications 
are also the tool to find the maximum period and relaxation. The results are summarized 
in Table [ITI] with those obtained from the eigenvalue analysis. The suppressed maximum 
period, which means the maximum period obtained by the direct matrix multiplications, 
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are just a half of those by the eigenvalue analysis. On the relaxation, the same situations 
happen except for the n N = 1 cases. The suppressed maximum periods are plotted in Figs. |2| 
and [| 

As the results from the suppression of the maximum values of the period and relaxation, 
we can find a simple relation of the period of rule 90 CA with periodic and Dirichlet bound- 
aries. The period of iV-site rule 90 cylindrical CA is half of N — 1-site rule 90 CA with 
Dirichlet boundaries. No such simple relations are found for rule 150 cases. 



IV. SUMMARY AND DISCUSSION 

The method of the eigenvalue analysis was applied to linear cellular automata with 
periodic boundaries. The trajectories of the states of CA were classified. The maximum 
values of the period and relaxation were evaluated by the eigenvalue equations. For rule 
90 simple relation was found between the eigenvalue polynomials with periodic 

boundaries and those with Dirichlet boundaries. 

The cylindrical CA enjoy the translational invariance with periodic boundaries. This 
translational symmetry suppresses the maximum values of the period and relaxation. 

In the previous papers [|12|,|13|] I had found the symmetry of the initial states suppress 



the period. So it is natural that the translational symmetry suppresses the period and 
relaxation. The translational symmetry changes the form of the transfer matrices U and 
consequently the form of the eigenvalue polynomials. The eigenvalue polynomials, however, 
do not tell us the suppression of the period and relaxation. 

The eigenvalue analysis may not be able to handle the translational invariance well. The 
matrix U for rule 90 CA can be decomposed as XJ = L + R, where L and R denote the 
translation operation to the left and right direction respectively. For rule 150 CA it can be 
describe as U = L + R + I. Namely the dynamics is described only with the composition of 
the translation operation. This may be the reason of the different results from the eigenvalue 
analysis and the direct matrix multiplications. 
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Let me investigate the cas of iV = 4 rule 90 CA with periodic boundary conditions for 
instance. The eigenvalue equation is A 4 = 0. This can be expressed with the transfer matrix 
U as 

U 4 = (L + R) 4 = L* + R 4 = I + I = 0, (4.1) 
where all operations are carried over GF(2). On the other hand, another identity 

U 2 = (L + Rf = L 2 + R 2 = (4.2) 
holds by virtue of the identity 

L m = R N- m ( 43 ) 
This identity may not be taken account into the eigenvalue analysis. 
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APPENDIX A: CLASSIFICATION OF ORBITS FOR DIRICHLET BOUNDARY 

CASES 

I had given the classification of orbits of rule 90 and rule 150 finite linear cellular automata 
with Dirichlet boundaries in paper II. Here I give brief proofs of the classification. 
For rule 90 cases, the eigenvalue polynomials obey the following recursion relation: 

dT\x) = xd^W - dSW- (Al) 

The equation ( |A1| ) is simplified for A = case to D^ OD \o) = /^^(O). With the special 
values D^ OD \o) = and Df OD \o) = 1, one can find simple relations: 
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DT\0) = D«(0) = M 9 ° D) (0) = ■•■ = !, (A2) 
D%$ (0) = « (0) = ^ 9 ° D) (0) = • • ■ = 0. (A3) 

The eigenvalue polynomials for even N have a constant term (Eq. (|A2|) ) . This means that 
rule 90 CA with Dirichlet boundaries with even number of sites show perfectly periodic 
motions. On the contrary, Eq. ( |A~3| ) shows that rule 90 CA with Dirichlet boundaries with 
odd number of sites have periodic orbits with relaxation paths or paths drawn into a null 
state. 

The another recursion relation 

< D) (A) = A 2 <_? (A) + DgS (A) ( A4) 

is derived from Eq. 0. By solving Eq. flA4]) with Lf 0D) (A) = A 3 and Z)f D) (A) = A 5 + A, 
a simple relation 

^ + D i ) (A)=ADf D )(A 2 ) (A5) 



is derived |15| . This gives explicit expressions of the eigenvalue polynomials for N = 2 n — 1 

cases: 

DS(A) = XD$£U#) =■= A 2 "" 1 - (A6) 
On the other hand, Eq. (|A1|) gives the explicit expressions of the eigenvalue polynomials 



with binomial coefficients as 

2 n-l 



«(A) = £ ( ^ - 1 - 3 j mod 2j A 2 " -1-23 '. (A7) 



Comparing Eqs. ( |A6|) and (|A~7l), identities on binomial coefficients 

mod 2 = (j ^ 0) (A8) 



2 n - 1 -j 
J 



are obtained. The eigenvalue polynomials for odd 2 n — 1) can be expressed with those 
for even A" by repeated usages of Eq. ( |A5| ) . So they can not be nilpotent. Therefore the 
following classification of orbits 
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U Un = I (N is even) , (A9) 

U Un+7TN = U 7TN (N is odd(^2 n -1)), (A10) 

U N = (N = 2 n -1), (All) 

is proven. 

For rule 150 cases, the recursion relation changes to 

Df 0D) (A) = (1 + \)D^f\\) - D%™\\). (A12) 

This leads to a simple relation D^ 50D) (0) = D ( ^\o) with D$ 150D) (0) = I^ 150D) (0) = 1 and 
L>5 150D ^(0) = 0. Therefore one canobtain simple relations 

^n OD) (0) = --- = J Df OD) (0) = l, (A13) 
^ 5 + °? ) (0) = --- = M 150D) (0) = l, (A14) 
^ 5 + °?(0) = --- = ^ 150D) (0) = 0. (A15) 

Orbits of rule 150 finite linear CA are perfectly periodic except N = 3n + 2 cases. 

One can prove that the eigenvalue polynomials of rule 150 CA with N = 3n + 2 are not 
nilpotent for n > 0. Corresponding to Eqs. (|A4]) and (|A5|) , recursion relations 



<° D) (A) = (1 + A 2 )Dj?(A) + 4?(A), (A16) 
C(A) = (l + ^i 150D) (A 2 ), (A17) 

hold. Applying them to odd 3n + 2, namely odd n, one can obtain the new recursion relation 

D { ^\\) = (1 + A)< n °\ 2)+2 (A 2 ). (A18) 

By the special case D^ 50D \\) = A 5 + A 4 , D^^W can n °t be nilpotent for odd n. For the 
even 3n + 2 cases, the new recursion relation 

Dg 50D) (A) = A 4 Dr 4 D) (A) + (1 + A 2 )D<S D) (A) (A19) 

is used. The lowest order term of Z^ 150D) (A) = A 8 + A 6 + A 2 is A 2 . And if D^\\) has a 
A 2 as the lowest order, 
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^ 5 + ° 8 D) (A) = A 4 DS 4 D) (A) + (1 + A 2 )D£f + °? (A) (A20) 

tells that also the lowest order of D^ r ^\\) is A 2 . Namely D^^\\) must have A 3n+2 and 
A 2 terms at least and can not be nilpotent for even 3n + 2(n > 2). 

Let me investigate special cases with N = 2™ — 1. The recursion relation Eq. ( A12 ) gives 



of D) W = E ( L< "g /2J ( N J j ) ( N ~ k 2i ) mod 2 ) A *- (a-) 

For special cases N — 2" — 1, this reduces to 

JV 

Z#-T(A) = £A fc , (A22) 

A;=0 

with the identity Eq. (|£8j) and ( 2n fc _1 ) mod 2 = 1(2™ - 1 > k > 0). This gives explicit 

expressions of the periods as IT^-i = 2 n . Therefore the classification of orbits for rule 150 
CA with Dirichlet boundaries 

ljn N = I (jV^3n + 2), (A23) 

U N = 1 (N = 2 n — 1), (A24) 

jju N+ n N = u* N (N = 3n + 2), (A25) 

is proven for N > 2. 
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FIGURES 

FIG. 1. Schematic features of the trajectories of cellular automata: (a) simple periodic orbit, 
(b) periodic orbit with relaxation, (c) limit point. 

FIG. 2. Actual period for rule 90 cylindrical cellular automata 

FIG. 3. Actual period for rule 150 cylindrical cellular automata 
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TABLES 



TABLE I. Eigenvalue polynomials of rule 90 cylindrical cellular automata 



A 

4 


A 


•5 


A + A + A 


6 


A + A 


7 


A' + A° + A 


8 


\8 
A 


9 


A y + A' + A + A 1 


10 


\ 10 i \6 i \2 
A + A + A 


11 


\11 i \ 9 i \ R i \ 3 i \1 

A 11 + A a + A° + W + A 1 


12 


\ 12 i \4 

A + A 


13 


\13i \lli \9i \3i \1 

A + A + A + A + A 


14 


\ 14 i \ 10 i \2 

A + A + A 


15 


\15 i \ 13 i \9 i \1 

A + A + A + A 


lb 


\ 16 
A u 


17 


\17 i \15 i \ 13 , \ 9 , \ 1 

A + A + A + A + A 


18 


\ 18 i \ 14 i \ 10 i \2 

A + A + A + A 


19 


\ 19 i \ 17 i \ 13 i \11 i \9 i \3 i \1 

A a + A ' + A ' + A + A a + A J + A 


20 


\20 i \ 12 i \4 

A + A + A 


21 


\21 i \ 19 i \ 17 i \11 i \9 i \5 i \3 i \1 
A + A + A + A + A + A + A + A 


22 


\22 i \ 18 i \ 10 i \6 i \2 

A + A ° + A u + A D + A 


23 


A 23 + A 21 + A 17 + A 9 + A 7 + A 5 + A l 


24 


A 24 + A 8 


25 


A 25 + A 23 + A 21 + A 17 + A 7 + A 5 + A 1 


26 


A 26 + A 22 + A 18 + A 6 + A 2 


27 


A 27 + A 25 + A 21 + A 19 + A 17 + A 5 + A 3 + A 


28 


A 28 + A 20 + A 4 
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29 A 29 + A 27 + A 25 + A 19 + A 17 + A 3 + A 1 

30 A 30 + A 26 + A 18 + A 2 

31 A 31 + A 29 + A 25 + A 17 + A 1 

32 A 32 
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TABLE II. Eigenvalue polynomials of rule 150 cylindrical cellular automata 

3 A 3 + A 2 

4 A 4 + l 

5 A 5 + A 4 + A 3 + A 2 + A 1 + 1 

6 A 6 + A 4 

7 A 7 + A 6 + A 3 + A 2 + A 1 + 1 

8 A 8 + l 

9 A 9 + A 8 + A 7 + A 6 + A 3 + A 2 

10 A 10 + A 8 + A 6 + A 4 + A 2 + 1 

11 A 11 + A 10 + A 5 + A 4 + A 1 + l 

12 A 12 + A 8 

13 A 13 + A 12 + A 11 + A 10 + A 9 + A 8 + A 5 + A 4 + A 1 + 1 

14 A 14 + A 12 + A 6 + A 4 + A 2 + 1 

15 A 15 + A 14 + A 11 + A 10 + A 9 + A 8 + A 7 + A 6 + A 3 + A 2 

16 A 16 + 1 

17 A 17 + A 16 + A 15 + A 14 + A 11 + A 10 + A 9 + A 8 + A 7 + A 6 + A 3 + A 2 + A 1 + 1 

18 A 18 + A 16 + A 14 + A 12 + A 6 + A 4 

19 A 19 + A 18 + A 13 + A 12 + A 11 + A 10 + A 9 + A 8 + A 5 + A 4 + A 3 + A 2 + A 1 + 1 

20 A 20 + A 16 + A 12 + A 8 + A 4 + 1 

21 A 21 + A 20 + A 19 + A 18 + A 17 + A 16 + A 11 + A 10 + A 3 + A 2 

22 A 22 + A 20 + A 10 + A 8 + A 2 + 1 

23 A 23 + A 22 + A 19 + A 18 + A 17 + A 16 + A 9 + A 8 + A 1 + 1 

24 A 24 + A 16 

25 A 25 + A 24 + A 23 + A 22 + A 19 + A 18 + A 9 + A 8 + A 1 + 1 

26 A 26 + A 24 + A 22 + A 20 + A 18 + A 16 + A 10 + A 8 + A 2 + 1 

27 A 27 + A 26 + A 21 + A 20 + A 17 + A 16 + A 11 + A 10 + A 3 + A 2 

28 A 28 + A 24 + A 12 + A 8 + A 4 + 1 
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29 A 29 + A 28 + A 27 + A 26 + A 25 + A 24 + A 21 + A 20 + A 17 + A 16 + A 13 + A 12 

+A 11 + A 10 + A 9 + A 8 + A 5 + A 4 + A 3 + A 2 + A 1 + 1 

30 A 30 + A 28 + A 22 + A 20 + A 18 + A 16 + A 14 + A 12 + A 6 + A 4 

31 A 31 + A 30 + A 27 + A 26 + A 25 + A 24 + A 23 + A 22 

+A 19 + A 18 + A 15 + A 14 + A 11 + A 10 + A 9 + A 8 + A 7 + A 6 
+A 3 + A 2 + A 1 + 1 

32 A 32 + 1 
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TABLE III. Period for rule-90 and 150 


cylindrical CA. 






Rule-90 




Rule-150 




N 


Polynomials 


Actual 


Polynomials 


Actual 


4 


U 4 = 


U 2 = 


U 4 = I 


U 2 = 1 


5 


U 7 = U 


u 4 = u 


U 6 = I 


U 3 = I 


6 


u 6 = u 2 


u 3 = u 


u 6 = u 4 


U 3 = U 2 


7 


u 15 = u 


u 8 = u 


U 14 = I 


U 7 = I 


8 


U 8 = 


u 4 = o 


U 8 = I 


U 4 = I 


9 


u 15 = u 


u 8 = u 


u w = u 2 


u 8 = u 


10 


u 14 = u 2 


u 7 = u 


U 12 = I 


U % = I 


11 


u 63 = u 


u 32 = u 


U &2 = I 


U 31 = I 


12 


u 12 = U 4 


u 6 = u 2 


u 12 = u 8 


u 6 = u 4 


13 


u 127 = u 


u 64 = u 


U 42 = I 


U 21 =I 


14 


u 30 = u 2 


u 15 = u 


U 28 = I 


U 14 = I 


15 


u 31 = u 


u 16 = u 


u 32 = u 2 


u 16 = u 


16 


U w = 


U 8 = 


U 16 = I 


U 8 = I 


17 


u 31 = u 


u 16 = u 


U 30 = I 


U 15 = I 


18 


u 30 = u 2 


u 15 = u 


u 32 = u 4 


u 16 = u 2 


19 


u 1023 = u 


u 512 = u 


JJ1022 _ j 


U 5U = I 


20 


u 28 = u 4 


u 14 = u 2 


U 24 = I 


U 12 = I 


21 


u 127 = u 


u 64 = u 


u 128 = u 2 


u 64 = u 


22 


u 126 = u 2 


u 63 = u 


U 124 = I 


U 62 = I 


23 


u 4095 = u 


u 2048 = u 


jj4094 _ j 


U 2047 = I 


24 


u 24 = u 8 


u 12 = u 4 


u 24 = u w 


u 12 = u 8 


25 


u 2047 = u 


u 1024 = u 


jj2046 _ j 


jj1023 _ j 


26 


u 254 = u 2 


u 127 = u 


U 84 = I 


U 42 = I 


27 


u 1023 = u 


u 512 = u 


u 1024 = u 2 


u 512 = u 
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28 U m = U 4 U 30 = U 2 U 56 = I U 28 = I 

29 f/32767 _ jj ^16384 _ jj jj32766 _ j jyl6383 _ j 

30 U 62 = U 2 U 31 = U U 64 = U 4 U 32 = U 2 

31 U 63 = U U 32 = U U 62 = 1 U 31 = I 

32 U 32 = U w = U 32 = I U w = I 
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